Dissipative Dynamics of Quantum Vortices 
in Superconducting Arrays 
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We consider a two-dimensional array of ultra-small superconducting grains, weakly coupled by 
Josephson junctions with large charging energy. We start from an effective action based on a 
microscopic tunneling Hamiltonian, which includes quasiparticle degrees of freedom, and study the 
resulting dissipative dynamics of quantum vortices. The equation of motion for a single vortex is 
deduced, and compared with a commonly adopted phenomenological model. 
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Vortex dynamics plays an essential role in understand- 
ing the transport properties of superconducting systems 
in external magnetic fields. For instance, it is closely 
related to the Hall resistivity Q and low-temperature 
magnetic relaxation Q , which has drawn much attention 
in the properties of vortices in two-dimensional (2D) and 
highly anisotropic three-dimensional superconductors || . 
In addition, the dynamics of vortices has been studied ex- 
tensively in superconducting arrays as well pj , for which 
the recent advancement in fabrication techniques allows 
one to control the physical quantities determining vor- 
tex dynamics, such as the vortex potential, the effective 
vortex mass, and viscosity. The superconducting arrays, 
therefore, provide a convenient model system, on which 
various theoretical predictions can be compared with ex- 
perimental results. Further, they may also shed light 
on physics of the high-temperature ceramic superconduc- 
tors, particularly in the polycrystalline form, which be- 
have in many respects like random arrangements of weak 
links. 

When the dimensions of the superconducting grains 
and the capacitances involved are small, the associated 
charging energy is non-negligible, and quantum dynamics 
of the phase comes into play at the macroscopic level fjl . 
In such an array of ultra-small junctions, the vortices, 
which are defined on plaquettes of the lattice, should be 
taken for quantum mechanical objects. Here it is gen- 
erally accepted that a vortex on a superconducting ar- 
ray is a rather well-defined point-like object with a finite 
effective mass and feels frictional force as well as 

non-dissipative transverse force in its motion, although 
there has been long-standing controversy as to the actual 
determination of the latter On the other hand, a 

vortex is a macroscopic object by nature, which raises the 
question regarding how to macroscopically quantize the 
proposed classical equation of motion for vortices in the 
presence of frictional force. It is a commonly adopted 
recipes for quantum dynamics of vortices to assume a 
frictional force linearly proportional to the vortex ve- 



locity and then phenomenologically quantize the result- 
ing equation of motion according to the Caldeira-Leggett 
procedure However, it is not obvious that the fric- 
tion should depend linearly on the vortex velocity, and 
the possibility of nonlinear behavior may not be excluded 
in advance. Indeed recent numerical simulations of the 
dynamics of an array containing one single vortex ap- 
pears to suggest the friction to be a nonlinear function 
of the vortex velocity p0| . This makes it necessary to 
investigate the quantum dynamics of vortices based on a 
model closer to the first principle, and desirable to obtain 
the effective action for vortices without phenomenologi- 
cal presumption. 

This paper presents an attempt toward such a goal: 
We start from an effective action based on a microscopic 
model for a 2D array of Josephson junctions consisting 
of ultra-small grains. In particular, we consider the case 
that the charging energy is non-negligible but still smaller 
than the Josephson coupling energy, so that vortices are 
well-defined, and show how the tunneling of quasiparti- 
cles introduces dissipation to the system. Using the dual 
form of the effective action, which describes the system 
of dissipative quantum vortices, we obtain the semiclas- 
sical equation of motion for a single vortex. We show 
that the damping on the vortex is is in general nonlin- 
ear in the vortex velocity and nonlocal in time. How- 
ever, it turns out that the nonlinear contribution in most 
cases becomes negligibly small except at very short length 
scales, and the frictional force in practice can be consid- 
ered to be linear in the vortex velocity, thus recovering 
the commonly adopted phenomenological model. 

An array of Josephson junctions can be described by 
the microscopic tunneling Hamiltonian pTlf] 



H = Hi + h t, 

1 (id) 



ff c>y> 



(1) 



where Hi represents the microscopic Hamiltonian, e.g., 
the BCS reduced Hamiltonian, for the zth grain. The 
coupling between neighboring islands results from the 
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transfer of electrons through the insulating barrier, de- 
scribed by HT,ij, and from the Coulomb interaction 
Hc,ij- Hx,ij is characterized by the Josephson coupling 
energy Ej whereas Uc,ij is characterized by the junction 
capacitance C, self-capacitance Co and the charging en- 
ergy Ec = e 2 /2C. Integrating out the quasiparticle de- 
grees of freedom leads to a macroscopic model, where the 
ith grain is characterized by the number m of the super- 
conducting electrons (Cooper pairs) and the phase fa of 



its superconducting order parameter jjjfl . The resulting 
Ambegaokar-Eckern-Schon (AES) model gives the parti- 
tion function in the form of a functional integral 



E 

{Mr)} 



2tt 



V[fa\ e^[-{So + S D )} 



(2) 



where the Euclidean action is given by 



f 13 1 
So = J dr n i<Pi + 2K E niC v> ~ K ^2 cos (<fe ~ A v ) 

~ (id) 

~<t>ij(r) - faj(r') 



r r r 

S D = dr I dr' a ( T ~ T ') { 1 
Jo Jo I 



(3) 
(4) 



with faj = fa — faj. Here we have rescaled the (imagi- 
nary) time in units of 1/uj p , where uj p = \J 8EqEj /ft 2 is 
the junction plasma frequency, and temperature in units 
of hujp. Further, we have introduced K = \J Ej / 8Ec 
and the dimensionless capacitance matrix 

Cij = (Cq/C + 4)<5ij- — — ^i,j-x — $i,j+y ~ $i,j-y- 

The bond angle Aij is given by the line integral of 
the vector potential due to the external magnetic field: 
= (27r/$o) /• A • dl, so that the plaquette sum 
gives the flux per plaquette in units of the flux quan- 
tum ($o = 2Trftc/2e) or gauge-invariant (magnetic) frus- 
tration, ^2 p A.^ — 27r/j, where i denotes the position of 
the plaquette. The quasiparticle degrees of freedom are 
effectively included through the damping kernel a(r) in 
the dissipative part Sd given by Eq. (||). In case that the 
grains are ideal BCS superconductors with energy gap A 
and normal state resistance Rn, the damping kernel is 
given by 

a(r) = ^JT?(A|t|) 



in the low-temperature limit (T — ► 0), where K\ is the 
modified Bessel function and we rescaled the gap energy 
A and normal state resistance Rn by A/ftaj p — > A and 
Rn/Rq — » Rn, respectively. It should be stressed here 
that the damping term entirely originates from the in- 
tergrain quasiparticle tunneling and includes neither the 
effects of the Cooper pair decay into the pool of normal 
electrons nor those of Ohmic shunt between grains. 

A variety of properties of Josephson-junction systems 
have been successfully described by the effective action 
in Eqs. (||) and (|J). We thus employ the AES model as a 
good starting point for an effective model for the dissipa- 
tive dynamics of quantum vortices. For this purpose, it is 
convenient to use the dual transformation, which rewrites 
the model in terms of the vortex variables instead of the 
original charge (Cooper pair) variables |0-[l5|. In the 
presence of the quasiparticle dissipation, the resulting ef- 
fective model reads [frfl 




with 



QV 

°o — 



dr 



- /,)Gy(nJ - fj) 
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ij 



(6) 
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S"d= dr dr' V a(r -t')\1- cos 
Jo Jo I 



/(-r)-^(r') 



(7) 
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where i,j now denote the dual lattice sites, i.e., the 
tilde signs representing the dual lattice sites have been 
dropped for simplicity. The vortex charge (in units of 
$o) n>i and the vortex phase §\ (of the macroscopic vor- 
tex wave function) are conjugate to each other, Gij is the 
lattice Green's function, and 9i = — 2irJ2j Gijrij. 

We now investigate the dynamics of a single vortex, 
focusing on the nature of the dissipation on the vortex. 



Such a single vortex may actually be introduced by ad- 
justing weak external magnetic field at zero temperature. 
For a vortex at position r(r), the dissipative action takes 
the form 



o D — 



r/3 

dr / dr' a(r 
i> Jo 



T')W(r(r)-r(r')), (8) 



where the non-local self interaction 



W(r(r)-r(r'))= ]T 



<rr'> 



1 — COS 



(r)-^(r') 



can be evaluated to reveal the logarithmic dependence 



^(r(r)-r(r'))«-ln 



r(r)-r(r') 



(9) 



(10) 



aside from an irrelevant additive constant. The pro- 
cedures of action minimization and analytic continua- 
tion O] then lead to the semiclassical equation of mo- 



tion describing the real-time dynamics of a single vortex 
driven by an applied current / (in units of the Josephson 
critical current Ij) in the x-direction 



2tt 



l Kr - 2 J dt' a(t - t')V r W{r(t) - r(t')) = 2irKIu x z. 



(11) 



The damping kernel a(t) in the real-time domain may be 
obtained by the analytic continuation 



Z7T 



(12) 



Alternatively, a(t) can also be obtained directly through 
the use of the real-time formalism PJT^], which, for the 
actions in Eqs. (§,§), gives a(t) related to the quasipar- 
ticle tunneling current iq P (w) (in units of Ij) via [ p"2]|T8"| 

ia(u>) = KI qp (u>). 



(13) 



Equation (|llj) is the most general semiclassical equa- 
tion of motion at length scales larger than the lattice 
constant. At zero temperature (T — 0), the quasiparticle 
tunneling current is simply given by 



4pM 



0, \u>\ < 2A 
27a;, \Q\ > 2A 



(14) 



where 7 = 1/uj p RnC. This allows one to write the equa- 
tion of motion in a more explicit and appealing form. 
The damping kernel in this case reduces to 



a(t) = 2lK— \5{t) - Asinc(2Ai)] 



(15) 



where sincx = (2/ttx) sin 2. Here it should be noted that 
the short time behavior of Eq. ( |l5j ) is valid only approx- 
imately, because of the high frequency cutoff in a(u)). 
With this simple form of the damping kernel, we finally 
obtain 



r + 7r — 7- 



dt' sinc[2A(t - t')]W"{r(t) - r{t'))r{t') 



-If, 



which, at long time scales, Eq. (|l^) takes the more explicit form 

f + 7± - 7A / dt' sinc[2A(t - t') 



n\r(t)-r(t'W 



-i/y. 

7T 



(16) 



(17) 
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The semiclassical equation of motion given by Eq. Jl6| ) 
or Eq. ( jl7| ) possesses two damping terms: One is lin- 
ear, but the other is nonlinear in the vortex velocity and 
nonlocal (memory-dependent) in time. When the vortex 
velocity is large, we have r(t')/\r(t) — r(t')| 2 ~ l/v(t') 
in Eq. (|l7|), and the nonlinear term becomes sufficiently 
small compared with the ordinary linear friction term. 
For small velocities, on the other hand, one must be care- 
ful about the short-wavelength cutoff [jl9| present in the 
function W(r). (Note that the continuum approximation 
has been used in the derivation of the equation of mo- 
tion.) To examine the behavior at low velocities, we have 
thus numerically integrated the equation of motion, and 
display the obtained frictional force as a function of the 
vortex velocity in Fig. |l|. It is revealed that the frictional 
force slightly deviates from the linear behavior at veloci- 
ties smaller than v c = A. It is of interest to note that v c 
may be written in the form v c = auoj /tt (in natural unit) , 
which corresponds to vortex hopping by one lattice con- 
stant during the characteristic time 1/loj = Ti/2\c\RnIj 
associated with the Josephson oscillation in a resistively 
shunted Josephson junction. 

The above analysis demonstrates that the frictional 
force on a vortex is practically linear in the vortex veloc- 
ity, in particular in the long-time and long-wavelength 
scale, where the semiclassical equation of motion is 
mostly concerned. Neglecting the nonlinear friction and 
rescaling the time in units of l/wj, we have Eq. jl^ ) in 
the reduced form 

ir/3 c r + Trr = -If, (18) 

which precisely corresponds to the commonly adopted 
phenomenological equation of motion describing the re- 
sistively and capacitively shunted junction (RCSJ) model 
with the Stewart-McCumber parameter (3 C = lojRnC 
H . This is remarkable in view of the fact that we have 
considered only quasiparticle tunneling, and neither the 
Ohmic shunt nor any other local damping sources have 
been included. 

In conclusion, we have considered a microscopic model 
for a two-dimensional array of Josephson junctions, in- 
cluding the quasiparticle degrees of freedom. From the 
effective action, which has been obtained without any 
phenomenological presumptions, the semiclassical equa- 
tion of motion for a single vortex has been deduced. It 
has been revealed that the quasiparticle tunneling pro- 
duces friction on the vortex motion. It includes a non- 
linear term in addition to the ordinary linear term al- 
though the nonlinear friction is in most cases dominated 
by the latter. At finite temperatures, the quasiparticle 
tunneling current /q P (a>) is smoothed out and approaches 
the Ohmic behavior (i.e. Iq P (u>) cx uj). In consequence, 
the non-linear friction term would become even less pro- 
nounced and essentially negligible. Similarly, arrays 
of Josephson junctions between d-wave superconductors 
are expected to display essentially linear behavior since 
d-wave superconductors have nodes in the momentum 



space at which the energy gap vanishes. This leads to the 
quasiparticle tunneling current with no sharp threshold 
in lo, even at T = 0, separating the high frequency Ohmic 
behavior and the low frequency non-Ohmic behavior [B0| . 
It is of interest to compare the nonlinear behavior found 
in this work with that obtained in Ref. In the lat- 

ter, numerical integration of the phenomenological RCSJ 
model has led to effective damping in the vortex motion, 
which becomes nonlinear as the velocity increases. Here, 
unlike the existing phenomenological approach, we have 
started from a microscopic model, and derived explicitly 
the equation of motion for a vortex in the system. The 
resulting damping term displays nonlinearity mainly in 
the low-velocity regime, which is in contrast with that 
obtained in the phenomenological approach. In addition 
to the quasiparticle tunneling, spin-wave excitations may 
provide another mechanism for the damping of vortex 
motion pl|j2^] . Since the spin- wave excitation requires 
energy of order of Tuo p [H for fiuj p > A, it is rather 
irrelavant. Furthermore, in the case of discrete charge 
states considered here, the region where the spin-wave 
damping can be ignored gets wider S| . Nevertheless for 
Tiijp <C A, however, contributions of spin-wave excita- 
tions may not be disregarded, and the investigation of 
the dissipation due to both the quasiparticle and spin- 
wave excitations will be a challenging topic. Finally, we 
remark that we have considered only the single vortex 
motion. In the case of many vortices, the vortex-vortex 
interaction becomes crucial, especially for Co <C C, which 
yields the logarithmically long interaction range. The in- 
teraction effects, based on the highly suggestive effective 
action in Eqs. (||) and (||) , will be an interesting topic for 
further study. 
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FIG. 1. Behavior of the frictional force (in arbitrary units) 
as a function of the vortex velocity. The dotted line represents 
the usual linear frictional force. The logarithmic scale should 
be noticed. The values of the parameters are Rn = 10.0 and 
KA = 1/8. 
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